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We obtain the rich phase diagrams in the Hubbard model on the triangular Kagome lattice as a 
function of interaction, temperature and asymmetry, by combining the cellular dynamical mean-held 
theory with the continuous time quantum Monte Carlo method. The phase diagrams show the asym- 
metry separates the critical points in Mott transition of two sublattices on the triangular Kagome 
lattice and produces two novel phases called plaquette insulator with an obvious gap and a gapless 
Kondo metal. When the Coulomb interaction is stronger than the critical value U c , a short range 
paramagnetic insulating phase, which is a candidate for the short rang resonating valence-bond spin 
liquid, emerges before the ferrimagnetic order is formed independent of asymmetry. Furthermore, 
we discuss how to measure these phases in future experiments. 
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Geometrically frustrated systems have shown a lot of 
interesting phenomena, such as the spin liquid and spin 
ice [3-li]. The charge and magnetic order driven by in- 
teraction is still one of the central issues in the field 
of strongly correlation systems 0-[l2|- Recently, a new 
class of two-dimensional materials CugX2(cpa)e ■ XH2O 
(cpa=2-carboxypentonic acid, a derivative of ascorbic 
acid; X=P, CI, Br) has been found fl3T4l5l ]. which is 
formed by an extra set of triangles (B-sites in Fig. 1) 
insides of the Kagome triangles (A-sites in Fig. 1). The 
Cu spins form a new type of geometrically frustrated lat- 
tice called triangular Kagome lattice (TKL). This lattice 
can also be realized by cold atoms in the optical lattices, 
in which the interaction between trapped atoms can be 
tuned by Fcshbach resonance and the kinetic energy can 
be adjusted by the lattice depth [l6l-22|. Although the 
effective spin models on this system have been investi- 
gated @, [23j , the real charge dynamics with spins and 
the phase diagram on this lattice have not been stud- 
ied. Different with other frustrated system, such as the 
triangular lattice and Kagome lattice, the "triangles-in- 
triangles" structure on the TKL induces two different 
sublattice. Therefore, it is desirable to investigate the 
Mott and magnetic transition on the TKL under the in- 
fluence of asymmetry which is induced by different hop- 
pings between two sublattices. 

Recently, many analytical and numerical methods have 
been dev elop ed to investigate the strongly correlated sys- 



tem [2443 ll. such as the dynamical mean-field theory 



(DMFT) [32|. However, the DMFT works ineffectively 
in the frustrated systems because the nonlocal correla- 
tion and spatial fluctuations are ignored. Therefore, a 
new method called cellular dynamical mean-field theory 
(CDMFT) has been developed to incorporate spatially 
extended correlations and geometrical frustration in the 
framework of DMFT [33rl36| by mapping the original lat- 
tice problem onto an effective cluster model coupled to an 



effective medium. The continuous time quantum Monte 
Carlo method (CTQMC) [37J is employed as an impu- 
rity solver in the CDMFT loop, which is more accurate 
than the translational QMC method due to the absence 
of Trotter decomposition. So, the CDMFT combining 
with the CTQMC can provide useful numerous insights 
into the phase transitions in the frustrated system [s , 34 [ . 

In the present Letter, we employ the CDMFT com- 
bining with the CTQMC to investigate the influence of 
asymmetry on the Mott and magnetic phase transition 
on the TKL based on Hubbard model. Two novel phases 
called plaquette insulator and Kondo metal induced by 
asymmetry is found in the phase diagram and their prop- 
erties is studied by determining the momentum resolved 
spectrums. By defining a magnetic order parameter, we 
can characterize a ferrimagnetic order for the large inter- 
action on the TKL. The phase diagram with the competi- 
tion between interaction and asymmetry is provided. In 
addition, the spectrum functions on the Fermi surface for 
the different interaction and asymmetry are presented. 
These interesting phases can be probed by the angle- 
resolved photoemission spectroscopy (ARPES) (38j, neu- 
tron scattering, nuclear magnetic resonance (NMR) and 
other experiments. 

We consider the standard Hubbard model on the TKL, 



H 



<ij>(7 



u ^2 n tt n ti + m 



(1) 



where tjj is the nearest-neighbor hopping energy, U is the 
Coulomb interaction, fj, is the chemical potential, cf a and 
Ci,j denote the creation and the annihilation operators, 
and riicr = c^ a Ci a corresponds to the density operator. 
The asymmetry factor is defined as A = t a b/tbb, which 
can be adjusted by suppressing the samples in experi- 
ments. For the convenience, we use tbb = 1.0 as the en- 
ergy unit. As shown in Fig. 1(a), there are two different 
sublattices on the TKL called A-sites and B-sites. The 
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FIG. 1: (Color online) (al) The unit cell of the triangular 
Kagome lattice (TKL) without asymmetry (A = 1.0). The 
open circles denote A-sites and the solid circles denote B- 
sites. The blue lines represent the hopping between A-sites 
and B-sites. And the red lines denote the hopping between 
B-sites. (a2) When A > 1.0, the TKL is similar with the 
Kagome lattice. (a3) When A < 1.0, the TKL is trans- 
formed into a system composed of many triangular plaque- 
ttes. (b) The thick red lines show the first Brillouin zone of 
the triangular Kagome lattice. The thin black lines corre- 
spond to the Fermi surface for the non-interacting case. The 
r, K, M, K' , M' points denote the points in first Brillouin zone 
with different symmetry, (c) The density of states of the tri- 
angular Kagome lattice for the A- and B-sites in the case of 
U = and A = t ab /t bb = 1.0. 



space group of the TKL is p6m same as the honeycomb 
lattice in the case of A = 1, as show in Fig. l(al). Fig. 
I(a2) shows the TKL is similar as the Kagome lattice 
when A > 1. Wc find the TKL can be transformed into 
a system composed of many triangular plaquettes in the 
case of A < 1, as shown in Fig. I(a3). The first Brillouin 
zone and Fermi surface of the symmetric TKL (A = 1) 
with U = is shown in Fig. 1(b). The density of states 
(DOS) for the different sublattices are shown in Fig. 1(c) 
in the case of U = and A = 1. They have the similar 
DOS around the Fermi surface. 

We improve the CDMFT by combining it with the 
CTQMC to investigate the correlation effects on the 
TKL. In the CDMFT, the original lattice is mapped onto 
an effective cluster model via a standard DMFT proce- 
dure. Thus, the effective medium g can be obtained by 
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FIG. 2: (Color online) Phase diagram of the triangular 
Kagome lattice at A = 0.6 with t bb — 1.0 as the unit of energy. 
The black solid lines with square points show the transition 
line of the A-sites, and the red dashed lines with circle points 
show the transition line of the B-sites. Two kinds of coexist- 
ing phases between red lines and black lines are the plaquette 
insulating phase and Kondo metal phase. Inset: The phase 
diagram of the symmetric TKL, i.e. A = 1, in which there are 
no coexisting phases. 



a Dyson equation: 

rV) = (E 



iu + [i — t(k) — Y,(iuj) 



)~ l +t{iw), (2) 



where t(k) is the hopping matrix element in the cluster, 
k is the wave vector within the reduced Brioullion zone 
based on clusters, S(iw) is the self-energy matrix, and uj 
is the Matsubara frequency. Then, we introduce an impu- 
rity solver, such as the CTQMC, to calculate the cluster 
Green function G(iuj). The new self-energy £(iw) can be 
calculated by Dyson equation E(iw) = g~ 1 (iuj)~G~ 1 (iuj) 
to close the self-consistent iterative loop. This CDMFT 
loop is repeated until the T,(iui) converges to a desired 
accuracy. The CDMFT contains the important informa- 
tion of nonlocal correlation and geometrical frustration, 
which is absent in the general dynamical mean-field the- 
ory (DMFT). In addition, the absence of Trotter decom- 
position in the CTQMC makes this method more accu- 
rate than the translational auxiliary- field QMC. 

The phase diagram obtained at A = 0.6 shows two 
coexisting phases and a reentrant behavior of the Mott 
transition when half-filling. The asymmetry can cause 
the separation of the phase transition points of the A- 
and B-sites. As shown in Fig. 2, with an increasing in- 
teraction at T = 0.5, the A-sites are found to translate 
from a metallic phase into an insulating phase and the 
B-sites stay in the insulating phase. This phase is called 
the plaquette insulating phase shown in Fig. 2. In this 
plaquette insulating phase, electrons are localized on the 
A-sites and the absence of next nearest neighbor hop- 
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FIG. 3: (Color online) The momentum resolved one-particle 
spectrum Afc(w) at A = 0.6 in the units of tbb = 1- (a) The 
metallic phase in the weak interaction case, such as U — 6.0 
and T — 0.5. (b) The Mott insulating phase in the strong 
interaction case, such as U — 9.0 and T — 0.5. An obvious 
single particle gap shows up around the Fermi energy, (c) 
The plaquatte insulating phase in the intermediate interac- 
tion case, in which the A-sites are insulating but B-sites are 
metallic, such as U = 7.6 and T = 0.5. A small gap shows up 
in this phase, (d) The Kondo metallic phase at U = 7.0 and 
T = 0.2, in which A-sites are metallic, B-sites are insulating. 
The single particle gap vanishes in this phase. 



pings causes the electrons can only be itinerant within 
the B-sites. When T < 0.34, the B-sites translates into 
the insulating phase, however the A-sites remain in the 
metallic phase, sec Fig. 2. This coexisting phase corre- 
sponds to a Kondo metal. The localized electrons on the 
B-sites act as the magnetic impurity, and the electrons 
on the A-sites are still highly itinerant. In this phase, 
the system shows the strong Kondo effect because of the 
high density of magnetic impurities. When the U > U c , 
for example, U c = 7.6 at T = 0.25, both the A-sites and 
B-sites translate into insulators. In Fig. 2, we find a 
reentrant behavior in the Mott transition of the A-sites 
caused by the frustration and the asymmetry, which is 
also found in the anisotropic triangular lattice [To| . This 
reentrant behavior divides the coexisting phases into two 
parts: the plaqucttc insulating phase and Kondo metal- 
lic phase. The inset figure in Fig. 2 shows the phase 
diagram of the TKL at A = 1.0. There are no coexisting 
phase and reentrant behavior found on the TKL when 
the asymmetry is absent, i.e. A = 1. 

In order to find the properties of the different phases 
in Fig. 2, we investigate the momentum resolved one- 
particle spectrum Ak(w) by using the Maximum Entropy 
Method (MEM) Hi). The results are shown in Fig. 3. 
The r, K, M and M' points have been defined in Fig. 1 
(b). In the case of T = 0.5 and U — 6.0, we find two 
obvious quasi-particlcs near uj = 4.0 and uj = 2.0, see 
Fig. 3 (a). The gapless behavior shows that the system 
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FIG. 4: (Color online) The evolution of double occupancy D 
on the A-sites as a function of interaction U with different 
temperatures at A = 0.6 in the units of tw, = 1. The inset 
figure shows the evolution of double occupancy on the B-sites. 
The arrows with different colors show the phase transition 
points at different temperatures. 



stays in the metallic phase shown in Fig. 2. When the 
interaction U increases, we find an obvious gap near the 
Fermi energy (see Fig. 3 (b)), which indicates the sys- 
tem becomes an insulator denoted in Fig. 2. Ak(uj) of 
the plaquette insulating phase is shown in Fig. 3 (c). 
In Fig. 3 (d), we find there is no obvious gap near the 
Fermi energy when the system stays at the Kondo metal 
state in Fig. 2. The momentum resolved one-particle 
spectrum function can be obtained by the ARPES ex- 
periment, and the plaquette insulator and Kondo metal 
phases can be found in certain samples. We define the 
double occupancy D on the TKL as 



D = dF/dU=±Y,( n ^l) 



(3) 



The decreasing of D on both sublattices shows that the 
suppressing of the itinerancy of electrons, which is a char- 
acteristic behavior of the Mott transition. When tem- 
perature drops, D increases due to the enhancing of the 
spin fluctuation at low temperature. The arrows indicate 
the phase transition point at different temperatures. The 
smoothly decreasing D characterizes a second order Mott 
transition when the interaction U increases. 

In Fig. 5, we show the evolution of spectrum function 
on Fermi surface as a function of k for different inter- 
actions U, temperature T and asymmetry A, which is 
defined as 



A(k;uj = 0) w - lim ImG k (uj + i0)/n, (4) 



where 



= Jf E e ik - (r -- 5) [w + M - *(k) - E(w)#. (5) 

and k is the wave vector in the original Brillouin zone. A 
linear extrapolation of the first two Matsubara frequen- 
cies is used to estimate the self-energy to zero frequency 
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FIG. 5: (Color online) The evolution of the spectrum function 
on Fermi surface in the units of tbb = 1. (a) A = 0.6, (b) 
A = 1.0, (c) A = 1.25. 

[3(| . In the case of A = 1.0, the spectrum function on 
Fermi surface shows six points near the M point in Fig. 
1 (b), similar as the non- interacting case. When inter- 
action U increases, the Fermi surface shrinks because of 
the localization of electrons. Fig. 5 (al) shows the Fermi 
surface is similar as the system composed of many trian- 
gular plaquettes at A = 0.6, due to the suppression of the 
hopping between the A-sites and the B-sites. The Fermi 
surfaces at A = 1.0 and A = 1.25 are similar as a Kagome 
lattice. When the interaction U keeps increasing, the 
Fermi surface is developing into a flat plane showing the 
localization of electrons, see Fig. 5 (a2), (b2), and (c2). 
The spectrum function can be measured by the ARPES 
experiment. 

In order to investigate the evolution of magnetic order 
on the TKL, we define a ferrimagnctic order parameter 
by: 

m = — ^2sign(i)({nit) - (n^)), (6) 

i 

where sign(i) = +1, for i=0, 2, 6; and sign(i) = —1, for 
i=l, 3, 4, 5, 7, 8. Fig. 6 shows the evolution of the single 
particle gap AE and m as a function of interaction U in 
the case of A = 1.0, T = 0.2. The Mott transition points 
of the A-sites and B-sites coexist when the asymmetry is 
absent. A gap opens in the case of U = 8.5. The ferri- 
magnctic order parameter m = at U < U c = 13.8 indi- 
cates the system is in a paramagnetic insulating phase. 
We argue that this paramagnetic insulating phase is a 
candidate for the short range RVB spin liquid due to the 
absence of long range correlations. When U > U c , the 
finite m means the system enters a ferrimagnetic state. 
The inset picture of Fig. 6 shows the evolution of single 
particle gap AE at A = 0.6 and T = 0.5, in which a gap 
opens on the A-sites first. 

Finally, we provide a phase diagram about the asym- 
metry A and interaction U at T = 0.2, see Fig. 7. 
As the asymmetry factor A increases, a coexisting zone 
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FIG. 6: (Color online) The evolution of single particle gap 
AE and ferrimagnetic order parameter m at A = 1.0 and 
T = 0.2 in the units of tn = 1.0. A paramagnetic metallic 
phase is found when the interaction is weak with AE = 
and m = 0. As the interaction U increases, a gap is opened 
and no magnetic order is formed with AE ^ and m = 0. 
This paramagnetic insulating phase can be a short range RVB 
spin liquid. An obvious magnetic order is formed when the 
interaction is strong enough with AE and m ^ 0. The 
insert picture shows the evolution of AE at A = 0.6 and 
T = 0.5. A plaquette insulator is found when the A-sites are 
insulating and the B-sites are metallic. 

containing the plaquette insulator and the Kondo metal 
shows up. This zone is suppressed from A = 0.9 to 
A = 1.11. When U > U c , such as U c = 13.8 in the 
case of A = 1.0, the system becomes a ferrimagnetic in- 
sulator with m ^= 0. Before entering the ferrimagnetic 
phase, a paramagnetic insulator state is found. We ar- 
gue this paramagnetic insulator state is a candidate of a 
short range RVB spin liquid due to the absence of any 
magnetic order and long range correlations. 

In summary, we have obtained the rich phase diagrams 
as the function of interaction U, temperature T, and 
asymmetry factor A. The asymmetry introduced in this 
strongly frustrated system can change the shape of the 
metal-insulator translation line and induces several in- 
teresting phases, such as plaquette insulator and Kondo 
metal. In the plaquette insulator phase, an obvious gap is 
found near the Fermi energy by investigating the momen- 
tum resolved spectrum function. The Kondo metal phase 
is an interesting phase where the partial sites work as the 
magnetic impurities in the metallic environment. We also 
find an interesting paramagnetic insulating phase which 
is related to the short range RVB spin liquid state before 
the ferrimagnetic order is formed. And the increasing of 
asymmetry does not suppress the emerging of this ferri- 
magnetic order. In addition, we have presented all kind 
of spectrum which can be used to detect these phases 
in real materials, such as CugX2(cpa)e • xH-^O, under 
applied pressure in experiments. Our studies provide a 
helpful step for understanding the coexisting behavior in 
metal-insulator transition, the formation of magnetic or- 
der and the emerging of spin liquid in frustrated systems 
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FIG. 7: (Color online) The phase diagram of the triangular 
Kagome lattice represents the competition between interac- 
tion and asymmetry for T = 0.2 in the units of tbb = 1. The 
region between the black lines with square points and the red 
lines with circle points denotes the coexisting zone which con- 
tains the plaquette insulator and the Kondo metal parts. A 
wide paramagnetic insulating region is found with an inter- 
mediate U. The blue lines with triangular points show the 
transition point to the ferrimagnetic insulator with a clear 
magnetic order. Insert: (a) Dimers formed in the paramag- 
netic insulator, which is a candidate for the short range RVB 
spin liquid, (b) Spin configuration of ferrimagnetic insulator. 



with asymmetry. 
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